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We discuss possible sources of systematic errors in the computation of critical exponents by renormalization- 
group methods, extrapolations from exact enumerations and Monte Carlo simulations. A careful Monte Carlo 
determination of the susceptibility exponent 7 for three-dimensional self-avoiding walks has been used to test the 
claimed accuracy of the various methods. 



1. Possible sources of systematic errors in 
the computation of critical exponents 

The major problem in a high-precision deter- 
mination ol critical parameters is the presence 
of corrections to scaling. Indeed they are a pri- 
ori unknown and cause a systematic error which 
is usually difficult to detect. The problem ap- 
pears in extrapolations from exact enumerations, 
where one is forced to use higher-order differen- 
tial approximants which however give stable re- 
sults only for long series. The problem arises in 
Monte Carlo (MC) simulations as well. Because 
the correlation length is necessarily much smaller 
than the lattice size L to avoid finite-size effects, 
the simulation is done in a region of temperatures 
so far from the critical point that neglecting the 
unknown corrections to scaling can deeply affect 
the final estimates. A clear example of this ef- 
fect can be found in the history of the estimates 
of the critical exponent v for three-dimensional 
SAWs: the first MC simulations jij with quite 
short walks suggested an exponent in agreement 
with the Flory theory v = 3/5, while a recent 
simulation ^ with very long walks (correspond- 
ing to a correlation length of £ (/3) rj 340) gave the 
much lower estimate v — 0.5877 ± 0.0006 which 
is in good agreement with renormalization-group 
(RG) estimates. It is clear that only a careful 
analysis of the Monte Carlo data can prevent from 



underestimating the error bars. 

On the other hand, a possible source of sys- 
tematic errors in the field-theoretic approach to 
critical phenomena could be the presence of non- 
analyticities on the real axis in the /3-function. 
Consider the critical behaviour of the coupling 
constant and of the mass of a generic system. 
If one includes, beside the dominant correction 
to scaling with exponent Ai, generic subleading 
ones with exponents 2Ai, A2 and analytic cor- 
rections, etc., one finds for the /3-function 
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One can easily see that confluent singularities on 
the real axis arise naturally if no specific hypoth- 
esis on the form of the corrections to scaling is 
made. In this situation the usual summation 
method || based on a Borel transformation and 
a complex-plane mapping, in which the request 
of analyticity of the /3-function on the real axis 
plays a crucial role, will converge very slowly. 

Nickel Q] suggested a new type of analysis, con- 
sisting in fitting the /3-function with functions 
with the expected cut singularity on the real axis, 
such as hypergeometric functions. With this type 
of analysis one finds ||, for the 0(n) cr-model an- 
alytically continued to n = 0, that the zero of the 
/3-function is at g* m 1.39, which is a value sensi- 
bly lower than the one obtained using the usual 
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type of analysis, g* = 1.421 ± 0.008. This fact 
reflects itself on the estimate of the critical expo- 
nent 7 in the following way: 

1 . the well-known value given by the standard 
analysis j| is 

7 = 1.1616 ± 0.11 (g* - 1.421) ± 0.0004 (2) 

2. the one obtained with the analysis sug- 
gested by Nickel || is 

7 = 1.1569 ± 0.10 (g* - 1.39) ± 0.0004 (3) 

This last value is in perfect agreement with our 
MC simulation Q , as shown below. 

2. The Monte Carlo simulation 

One can hope to be able of discriminating 
between the two proposed types of analysis by 
means of a Monte Carlo simulation. The SAW 
is the best possible test-case because the simu- 
lations are not affected by finite-size effects and 
because there are algorithms which have autocor- 
relation times much smaller than those of the al- 
gorithms which are available for other systems. 
To give an idea, for our algorithm the autocorre- 
lation time in CPU units scales as r ~ £ 2 , which 
can be compared with the best algorithms for the 
Ising model, for which r ~ £, d+z , with z ps 0.4. 
A detailed description of the algorithm used in 
this work can be found in p] . The algorithm 
works in the ensemble of pairs of walks with fixed 
total length N to t- One can make inferences on 
the value of 7 from the observed distribution of 
length of one of the two walks, which is in our 
case 
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where cj^ 1 denotes the number of walks of length 

Nx and Z = J2nT=o c Ni c N tot -N 1 - 

Instead of making an a priori hypothesis on the 
form of the corrections to scaling to fit the MC 
data, we prefer to keep only the leading term in 
the scaling law for cn 



and make a type of analysis which is sensible to 
the presence of corrections to scaling, in the fol- 
lowing way. We progressively reject those walks 
with length shorter than a certain length N m i n 
and greater than N to t — N m i n . Of course the 
estimated values of 7, *f(Ntot, N m i n ), depending 
on the range of allowed lengths and on N to t, con- 
verge for N t ot — > co and N m i n — > 00 to the correct 
value. 

High-statistics runs at different total length, 
giving an idea of how strong the effect of cor- 
rection to scaling can be, are shown in Table [|. 
The run at N to t = 200 shows a strong depen- 
dence on the cut N m i n , but is in agreement, for 
large values of the cut, with |3j and with previ- 
ous MC simulations g. For N tot = 2000, the 
estimated values are lower and much flatter, but 
still sensibly biased by corrections to scaling. For 
the final estimate we used a weighted average 
of the estimates obtained using walks of length 
N tot = 20000 and N tot = 40000. The effective 
exponents are shown in Figure |l|, where the last 
point on the right denotes the final estimate with 
its error bar. The estimates reported in Fig. |l| 
show a very small dependence on N m i n but, of 
course, one cannot conclude that the estimates 
are not biased by corrections-to-scaling effects. 
Indeed to understand the residual systematic ef- 
fects we must compare estimates with different 
values of N to t- Assuming that the corrections to 
Eq. (||) scale with a subleading exponent A > 0.5, 
one confirms that our data at the larger values of 
Ntot have a systematic bias which is less than the 
statistical error. A very conservative way of de- 
termining the residual systematic error consists 
in making the hypothesis that the effective esti- 
mates 7 tend toward the true value following the 
law 
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Fitting our data to this law, requiring that (|6j) 
reproduces the highest estimated exponent for 
Ntot = 2000 and the lowest one for N to t = 40000 
with A = 0.5, we estimate B « 0.06. This result 
clearly leads to an overestimate of the systematic 
error if A > 0.5. With this type of analysis of the 
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Table 1 



Estimates of 


7 for Ntot = 200 and N tot = 2000. 
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the number of iterations. 
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Figure 1. Estimates of 7 obtained from the 
weighted average of the results with Ntot = 20000 
(number of iterations Nuer — 10 8 ) and N to t = 
40000 (Nuer = 8.5 • 10 8 ). For clarity the sta- 
tistical error is reported only in two cases. The 
dotted line and the rightmost point indicated by 
a diamond is our final estimate with its error bar. 



confluent singularities of the /3-function on the 
real axis at g*. We hope that further analytic 
work can be done in this direction to clarify this 
point. Previous Monte Carlo and exact enumer- 
ation determinations appear to be incorrect as 
well, showing how important and difficult is the 
determination of the systematic error due to cor- 
rections to scaling. 
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data, we estimate 

7= 1.1575 ±0.0006 . (7) 

3. Conclusions 

The present work indicates that the previous 
estimate of the critical exponent 7 obtained from 
the RG ^-expansion is significantly biased up- 
ward, possibly because of the presence of strong 



